x=-2*pi:2*pi/360:2*pi;
y=-2*pi:2*pi/360:2*pi;
[X,Y]=meshgrid(x,y);
Z=cos(X).*sin(Y)./Y;
plot3(X,Y,Z);
shading flat